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abstract 

Using a recently proposed new renormalization group method (tensor renormalization group), 
we analyze the Ising model on the 2-dimensional square lattice. For the lowest order approxi- 
mation with two domain wall states, it realizes the idea of coarse graining of domain walls. We 
write down explicit analytic renormalization transformation and prove that the picture of the 
coarse graining of the physical domain walls does hold for all physical renormalization group 
flows. We solve it to get the fixed point structure and obtain the critical exponents and the 
critical temperature. These results are very near to the exact values. We also briefly report the 
improvement using four domain wall states. 



1 Introduction 

Recently a new type of renormalization group method called the tensor renormalization 
group (TRG) is introduced by Levin and Nave [2] and has been applied to various models, 
particularly classical and quantum spin models in 2- dimension [3] S]. In this article we 
concentrate on the square lattice Ising model. The tensor renormalization group method 
for spin systems can be interpreted as the domain wall renormalization group. We hope 
it might give deep insight for the idea of the coarse graining of such topological objects 
like conserved walls. 

Our aim here is to clarify in detail the properties of this domain wall renormalization 
group transformation in the lowest order approximation with only two domain wall states 
on the coarse grained links. We will be able to write down the explicit analytic form 
of the renormalization group transformation (RGT) with the help of the physical region 
condition. This analytic formula will help us to prove that the picture of the corase 
grained domain walls does hold for all renormalization group flows starting from the 
physical points. 

Then we analytically prove the existence of a single non-trivial fixed point controlling 
the phase transition and evaluate the correlation length critical exponent which is actually 
unexpectedly good value, that is, within 2% error of the exact value. 

We also obtain the critical temperature within 10% error of the exact value. Us- 
ing this lowest order RGT we also numerically calculate the partition function for finite 
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Figure 1: Domain wall representation of the square lattice Ising model 

size systems as a function of the temperature. In the infinite size limit, we observe the 
logarithmic divergence of the specific heat in the neighborhood of the criticality. 

We finally mention quickly that enlarging the coarse grained state space we set up 
4-state version of the domain wall renormalization group and it gives an yet better result 
for the critical temperature within 2% error. 

We would like to stress here about new arguments and results in our study. We obtain 
explicit analytic forms of the renormalization group transformation and solve its fixed 
point condition to have a single non-trivial fixed point solution. Although we work in the 
lowest order approximation, this is a good example of partially solvable renormalization 
group equation to clarify total phase structure. Analytic forms give us much greater 
inspiration than numerical calculations in the black box, and we hope it will contribute 
to physical formulation of renormalization group equation for topological objects in the 
higher order approximation. In this course of calculation it is important that we define 
the physical region condition of parameters, because the condition assures that the notion 
of the conserved domain walls hold also for coarse grained variables. Due to this fact, the 
coarse grained variables can be regarded as to describe the coarse grained domain wall 
configurations, and therefore we can call the renormalization group as the domain wall 
renormalization group. Also the critical behaviors we evaluated for this system of the 
square lattice Ising model has not been reported explicitly elsewhere. 



2 Domain wall representation and RGT 

It is well-known that in the 2-dimensional Ising model, any spin configuration is equiv- 
alently represented by a domain wall configuration defined on the dual link of the dual 
lattice. As for the square lattice, the dual lattice is also square. Domain walls are said to 
exist on a dual link when two spins separated by the dual link are different (Fig{TJ). 
Note that a domain wall configuration represents two spin configurations of Z2 pair. Thus 
summing up all domain wall configuration gives just a half of the partition function and 
this factor 1/2 does not affect physical quantities discussed below. 
We define the partition function of the Ising spin system by 

Z = J2 ex P I PJ^Wj ) ■ (!) 

a \ n.n. / 

According to the recently proposed general method of the tensor renormalization group 
method[2], we express the partition function in terms of a product of vertices (sometimes 
called tensors, the origin of the name, tensor renormalization group). 

% — X! T abcd T ed f g ■ ■ ■ (2) 
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Figure 2: Integrate out S- vertex to define renormalized T- vertex 



The vertex T is defined on each dual site and has 4 indices representing domain wall state 
on 4 dual links. We assign for no-domain wall and 1 for existence of domain wall. Then 
the components of T are given by 

Toooo = exp(2/3J) , 

^oioi = ^oiio — 1 j (3) 
Tim = exp(-2/3J) , 

where the index cyclic symmetry is assumed. 

The coarse graining procedure consists of two steps. The first step is to break T- vertex 
into a product of two S"-vertex. From the Feynman diagram view, this step can be seen 
as introducing an intermediating boson for the 4-fermi weak interactions. 




The intermediating line will become a coarse grained dual link and coarse grained domain 
wall states are defined on it. If we take 4 states for that, we do not lose any information. 
However we like to set up a finite dimensional system of the renormalization group trans- 
formation and we have to discard some degrees of freedom here. Then we take only 2 states 
on the intermediate states. It means this decomposition of T-vertex loses information, 
and this is the approximation of this type of renormalization group transformation. 

The 2-state approximation must be the lowest approximation and we will analyze its 
RGT in detail. Which two states among totally 2x2 = 4 states should be picked up 
? There are two ways to determine it. One is a physical condition that coarse grained 
states should obey notions of the coarse grained domain walls. The other is a practical 
condition that we should minimize the information loss of this procedure. Fortunately, we 
prove below that these two conditions are consistent for some region of parameters which 
we call physical region, and it is respected by the physical finite temperature system. 

The second step of RGT is to integrate out all ^-vertices and will make the system 
described only by new variables, the intermediate states. The new variables are nothing 
but the coarse grained domain walls on the newly made coarse grained dual links (FigfSJ). 
These two steps define RGT with scale factor of \/2. 
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3 Analytic expression of RGT 



We express T-vertex (with four legs) by a 4 x 4 matrix. We combine two legs each into 
row or column index of the matrix. The correspondence rules are 



{00} =► {1} , {11} =► {2} , {01} =► {3} , {10} =► {4} 
Then we assume the texture of T matrix as follows: 
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(6) 



where a,b,c are real positive parameters obeying the conditions defined by 

< a < b < 1 , 
< c < b , 
2b < c+ 1 , 



(7) 



which we call the physical region condition. In general, renormalization group transforma- 
tion modifies T u component and we always normalize the matrix by dividing all elements 
by a constant so that Tu should be 1. 

These two properties, the texture and physical region condition, are satisfied by the 
initial T matrix given by 



a = b = a , 

2 



(8) 



exp(-2/3J) < 1. 



Now will prove that our RGT does not break these two properties and therefore all flows 
starting from the physical region have the same texture and belongs to the physical region. 
First we calculate eigenvalues of T. 



1 



^1,2 

A 3 a = a±b . 
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(9) 



We examine the absolute values of these 4 eigenvalues and select two larger components 
according to the policy of TRG. Using the physical region condition, we can prove the 
following inequalities 



Ai > b , A 3 > b , 
|A 4 | < b , |A 2 | < b 



(10) 



Thus we take two eigenvalues Ai, A 3 to define S vertex. Note that these two eigenvalues 
are positive. 
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Now we can split T-vertex into a product of two ^-vertices. First we introduce two 
orthogonal 2x2 matrices which diagonalize T as follows: 



T 



( 1 b \ 
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V b a ) 
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V A 4 J 



where R\ t 2 are defined by angles 9\^ 

Ri = 
R 2 = 



cos 9 1 — sin 9\ 
sin cos 6>i 

cos 6*2 — sin #2 
sin ^ 2 cos ^ 2 



Therefore we have the following decomposition: 

T = S A -S B , 

where S A is 4 x 2 matrix and S 113 is 2 x 4 matrix given by 



s 1 








Aici 




\ 


A3C2 

A3S2 / 
ATsi 

VA3C2 V^3 s 2 



R\ 

Rn 



11^ 



(12) 



(13) 



(14) 



It should be noted here that these S'-vertices satisfy S A = (5 B )* and they are essentially 
the same. This decomposition is nothing but the singular value decomposition which is 
frequently used in order to pick up important degrees of freedom. Hereafter we use simple 
notations as follows: 



si = sin 9\ , c\ = cos 9\ 



*2 



sin 9 2 , c 2 



cos 9o 



(15) 



4 Feynman rules to calculate RGT 

The T-vertex decomposition is expressed as 





(16) 



T S A S B 

where the inner thick line represents two possible states. We assign names to these two 
states to be (Ai component) and 1 (A3 component) respectively. This assignment is 
essential for interpreting the coarse grained domain wall variables, which will be clear 
soon. The ^-vertex gives the following Feynman rules: 
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where we draw single line for 0-state and double line for 1-state. These rule applies to 
all legs. Now we see S- vertex conserves the number of legs of double line (1-state), thus 
conserves the domain wall. 

The renormalization group transformation is defined by one-loop diagrams and they 
are calculated as follows: 
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(19) 



AiA 3 c 2 s 2 (c\ + s?) ,(20) 
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Note that the conservation of domain wall prohibits other diagrams. These amplitudes 
give the renormalized T-vertex. Taking account of the discrete rotational symmetry, 
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we have the renormalized T-vertex (T r ) 
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where we made additional total renormalization of the matrix elements so that T' X1 should 
be equal to 1. 

Thus we have proved that RGT conserves the texture defined in Eq.([6]). This comes 
from the property that the S-vertex conserves the domain wall with our assignment of 
coarse grained domain walls. Also the discrete rotational symmetry is necessary. In fact 
the physical region condition ([7]) is not necessary for this texture conservation. Strictly 
speaking the necessary and sufficient condition is that each of the two coarse grained 
components should come from the upper block and the lower block respectively. The 
physical region condition corresponds to another property that each of the two largest 
values among 4 eigenvalues are in the upper block and the lower block respectively. Thus 
the physical region condition assures that both of the best approximate coarse graining 
and the domain wall conservation holds simultaneously. 

We should discuss here the notion of coarse grained domain walls. The S'-vertex 
defines the relation of the micro domain walls and macro (coarse grained) domain wall. 
The micro domain wall does not always generates macro domain wall. The small loop 
of the micro domain wall does not correspond to the macro domain wall, whereas the 
domain walls connecting two coarse grained vertices (domain walls residing on the coarse 
grained links) are defined as macro domain walls. These interpretation naturally defines 
the coarse graining of the domain walls and renormalization of the T- vertex. 

5 Renormalization Transformation 

The renormalization transformation for parameters a, b, c are written down as follows: 



a 



2X 3 C1S1C2S2 

"a7 4 + 4 

A3 c 2 s 2 



(24) 



\i 4 + 4' 
2A3 c|s| 
1^4 + 4 ' 



Angles are obtained as a function of parameters a, b, c 
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Using these angles we rewrite the transformation 
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(26) 
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where eigenvalues and diagonalization angles are 



2b 

A 3 = a + b , s = . = . (27) 

'I -c) 2 + 4b 2 



This completes the explicit analytic formula of the domain wall renormalization group 
transformation. 

Finally we check the physical region condition for renormalized parameters. Due to 
the inequality < a < b, we have 

A 3 < 2b . (28) 
On the other hand, we define a function /(A) 

/(A) = A 2 -(l + c)A + c-6 2 = (29) 

where Ai is a larger root of this equation. We have the following inequality 

f(2b) = (1 - b)(c - b) + b(2b - c - 1) < (30) 

and it proves Ai > 2b. Then we have 

£<1. (31) 

Also using < s < 1, we can finally prove 

< a' < b' < 1 , 

< d < b , (32) 
2b' < c + 1 . 

Therefore we have proved that the physical region condition ((?]) is conserved by the 
renormalization transformation. Then the above analytic forms of the transformation 
and physical region condition always hold for any flows starting from a point in the 
physical region. 

6 Fixed points 

According to the standard procedure, we look for fixed points of the renormalization 
transformation. We can see there are two trivial fixed points. 

Low temperature fixed point:{a = 0, 6 = 0, c = 0} 

High temperature fixed point :{ a = 1, 6 = 1, c = 1} 
These fixed points are not in the physical region but on the boundary of the region. They 
are both infrared fixed points which almost all flows (except for those on the critical 
surface) in the physical region approaches in the infrared limit. We call them low or 
high temperature fixed point respectively because they are equal to the physical initial 
T- vertex with a = 0(/3 — > oo) or a = — > 0). 
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We check the eigenvalues of the renormalization transformation in the neighborhood 
of these trivial fixed points and find that all eigen values are less than unity for both 
fixed points. Thus, we have two infrared fixed points at the boarder of the physical 
region. Then there must necessarily exist a critical surface which divides the space into 
two subspace (different phases). Phases are characterized by these infrared fixed point. 
On the critical surface there must exist at least one (non-trivial) fixed point. 

The non-trivial fixed point satisfies the following set of equations 



A 



?> 



X x ) 2-s 2 ' 

Aj 2-s 2 ' 
; A^ 2 < 



(33) 



AW 2-s 
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where Ai,A3,s are defined in Eq. fl27j) . Noticing that parameters s and Ai are functions 
of b, c only, we first solve aD Dividing the first equation by the second equation in the 
above, we have 

2b 2 

a = bs= = . (34) 

^(1 - c) 2 + Ab 2 



Also dividing the second equation by the third equation in Eq. (j33[) . we have 

b 2 (2-s 2 ) = c. (35) 

and we solve it with respect to b, 



1 



6=-y2c-(l-c) 2 + ^/4c 2 + (l-c) 4 . (36) 

Next adding the first and second equations in Eq. fl33l we have 

A x (2-s 2 ) = 1 + s . (37) 

Substituting Eq.( l36|) into the above equation, we have an equation for c. 

136c 5 - 145c 4 + 116c 3 - 54c 2 + 12c - 1 = . (38) 

This 5th order polynomial equation is easily found to have only one real root. We finally 
get the unique non-trivial fixed point 

c* = 0.238902743 , 

b* = 0.402938077 , (39) 
a* = 0.292942734 . 



Numerical results for the critical surface and the renormalized trajectory are drawn 
in FigJSJ Note that we plotted here whole region of < a, b, c < 1. There is no singular 
behavior at the boundary of the physical region condition once the renormalization trans- 
formation is given. Also the physical initial points and a sample renormalization group 
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Figure 3: Phase structure of domain wall RGT 



flow are plotted there. The cross point of physical initial points and the critical surface 
determines the critical temperature. We have the criticality in terms of a 

a c = 0.37036 , (40) 

which should be compared with the exact value a c (Exact) = y/2 — 1 = 0.41421. Our 
result suffers about 10% error. 

We linearize the renormalization transformation around the non-trivial fixed point and 
find the eigenvalues 

{1.4224236,-0.29079698,0}. (41) 

There is only one eigenvalue A = 1.4224236 which is larger than unity. Thus the fixed 
point has only one relevant operator as is reduced by the general argument. According 
to the standard method [I] to evaluate the critical exponent v of the correlation length 
divergence, we have 

v = = 0.98357 . (42) 

log A 

Here we have used that the renormalization transformation changes the lattice scale by 
factor \/2. This result seems extremely good, less than 2% error compared to the exact 
value v = 1, considering that our domain wall renormalization group has only two states 
on the coarse grained link and only 3-dimensional interaction space has been used. 

We proceed to investigate the behavior of the mean energy and the specific heat. 
The partition function is calculated by tracing out the T-vertex. Note that we have 
to take account of the additional total normalization factor which is introduced at each 
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renormalization step to make T u to be unity. We denote the factor by C n which is used to 
divide T at the n-th renormalization step. After n-times renormalization transformation, 
the T-vertex is renormalized to become and we define 

a\ _ r< r r( n ) _ n 2 ™ r 2 * 1 ' 1 r< 2n ~ 2 C T^> (A "i\ 

Zj {p) — O tot l abab — O O x 2 • • • O n l abah , {46) 

which is the partition function of the system with size 2 n (the total number of sites). 

We calculate the mean energy and the specific heat per site. The results are shown 
in FigsJH and EJ In the mean energy plots, increasing n, there start to appear a would-be 
singular point near (3 J ~ 0.5. The specific heat would diverge at the point due to non- 
analyticity at n — > oo. The divergence behavior of the specific heat is almost completely 
logarithmic as is derived by the scaling relation. 

We can increase the number of states defined on the coarse grained link. As a next 
approximation we set up 4-state version of the domain wall renormalization group. The 
corresponding T-vertex has 4 4 = 256 components. The texture analysis shows that only 
25 components are non- vanishing in the procedure of the renormalization group transfor- 
mation. The details of our study will be reported in a full paper [5]. Here we only mention 
some results. The critical temperature is obtained as 

a c [4-state] = 0.42205 , (44) 

and the error is now less than 2%. This is a great improvement compared to the 2-state 
result. As for the correlation length critical exponent, we have 

^[4-state] = 0.98359 , (45) 

and this is almost equal to the 2-state result. 

In the 4-state version RGT, there are some subtle issues to be discussed in detail. For 
example, notion of the coarse grained domain wall is not trivial like 2-state version. To 
evaluate the magnetic quantities like the magnetization and the susceptibility, we need to 
introduce external source field. This will break the symmetry property we have used in 
this article and it needs additional consideration. 
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